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Factorization based on the soft-coUinear effective theory (SCET) can be used to reduce the number 
of hadronic parameters in an isospin analysis of B — > tttt decays by one. This gives a theoretically 
precise method for determining the CP violating phase 7 by fitting to the B nn data without 
C^o^o. SCET predicts that 7 lies close to the isospin bounds. With the current world averages we 
find 7 = 75° ± 2°^j3o , where the uncertainties are theoretical, then experimental. 



Measurements of CP violation are an important tool 
to look for physics beyond the standard model (SM) 0. 
Standard model measurements of CP violation in B- 
decays are usually expressed in terms of the angles a, 
fi, 7. To test the SM picture of CP violation one looks 
for inconsistencies by making measurements in as many 
decay channels as possible. 

Important observables for measuring 7 (or a) are 
the CP asymmetries and branching fractions in _B — s- 
TTTT decays. Unfortunately, hadronic uncertainties and 
"penguin-pollution" make the data difficult to interpret. 
Gronau and London (GL) '2*1 have shown that using 
isospin, Br(i3 tt+tt-), Br(B+ tt+i:°), By{B 
7r°7r°), and the CP asymmetries C^+.n._, Ct^Ot^q, 
one can eliminate the hadronic uncertainty and deter- 
mine 7. Thus data is used to determine the 5 hadronic 
isospin parameters. This year Babar and Belle Q re- 
ported a first observation of Ct^Qt^o. Unfortunately, the 
uncertainties in CttO^o and Br(i? — + tt^tt'') are still too 
large to give strong constraints, leaving a four-fold dis- 
crete ambiguity and a ±29° window of uncertainty in 7 
(at 1-cr) near the SM preferred value. 

In this letter we observe that the soft-collinear effec- 
tive theory (SCET) predicts that one hadronic pa- 
rameter vanishes at leading order in a power expan- 
sion in Aqcd/_B^, and that this provides a robust new 
method for determining 7 using the experimental value 
of /?. The parameter is e = Im(C/T), where T and 
C are defined below and are predominantly "tree" and 
"color suppressed tree" amplitudes. From we know 
that e vanishes to all orders in as{y/K^^Juy^) since the 
"jet-function" does not involve a strong phase, and so e 
receives corrections suppressed by AqcD/^'ir or as(mfc). 
Our method does not rely on a power expansion for any 
of the other isospin parameters. Thus, issues hke the 
size of charm penguins and whether "hard-scattering" 
or "soft" contributions dominate the B ir form fac- 
tors Is,, J., Jj 8.., ,9, ,KI, nj are irrelevant here. Our analy- 
sis remains robust if so-called "chirally enhanced" power 
corrections H are included. It differs from the QCDF 
and pQCD analyses; for example we work to all or- 
ders in Aqcd /iTT-b for most quantities and do not use QCD 
sum rule inputs. 
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For later convenience we define the ratios 
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0.446 ±0.064, 
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= 0.293 ±0.064. 
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and quote the product i?„ C^o^o = 0.082 ± 0.116. 

To obtain general expressions for these observables, we 
use isospin and unitarity of the CKM matrix to write 

A{B° ^n+n-) = e-'^ |A„| T - |A,| P 

+ (e-^^|A„|-|A,|)P,\,, 

A{B° ^ tt\'') ^ e~'^ \Xu\C + \Xc\P 

+ (e-^^|A„|-|A,|)(Fl-F,ij, 

V2AiB~ it\') = e-'-' (T + C) 



+ (e-''^|A„|-|A,|)P, 



2 

cw ■ 



(3) 



Here A„ = V^.V^^, Xc = V^.V;^. The CP conjugate am- 
plitudes are obtained from ^ with 7 —7. T, C, P and 
the electroweak penguin amplitudes Pqw^ complex. 

The amplitude is related to T and C by 

isospin [13 ■ An additional relation for P™ can be 
obtained using SCET at lowest order in A/E^^ and 
c^simb) 0. For the dominant coefficients Cg^io we find 



PpL = ei r + 62 c , p„^^ = 63 r + 64 c , 



(4) 



Clo(Ci — C3)+Cg(C4 — C2) 
(Ci + C2)(Ci-C2-C'3 + C4) ' 
Cg{Ci+C4) — Cio{C2 + C3) 
(Ci+C2)(Ci-C2-C3+C4) ' 

ea = e4 = ^(C9 + Cio)(Ci + C2)-i 
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FIG. 1: a) Isospin triangle in |Au| sector, and b) the rescaled 
triangle with solutions for positive and negative e shown. 



with Ci from the electrovi'eak Haniiltonian. Since 
e3|r|/|P| = e^(j>l + pI)^/^\K\/\K\ ^ 0.06 for typical 
values of the parameters Ps and pc (from below), we esti- 
mate that the electroweak penguins give at most a ~ 6% 
correction to any amplitude. It would be easy to include 
, but for simplicity we neglect them in what follows. 
SCET allows contributions from C7 and Cs to be in- 
cluded in igj), giving 63 = -1.5 x 10^^, 64 = -1.3 x 10"^. 

Of the five remaining isospin parameters, one, |A„(T + 
C)|, is fixed by Br(_B~ tt'^tt") and just sets the overall 
scale. We choose the remaining four parameters as 
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tr 



\M TP (P\ |Ac| , (P 

e= Imf — ^ . 
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The parameters pc and ps determine the size of the "pen- 
guin" contribution P relative to the "tree" T, and the 
parameters tc and e determine the shape of the isospin 
triangle as shown in Fig.Q] Relations to parameters used 
previously ^'^s = Pc+Ps tan (5c — Ps/pc- 

In terms of the parameters in lO the observables can 
be written as (neglecting electroweak penguins) 

S^+^- = -[sin(2/3+27) + 2sin(2/3-f7)pc 

+ sin(2/3) {pI +pI)\ [1 + 2pc cos -f+pj +pl] ' ] 
C - = sin 7 

l + 2pcC0S7-|-p2+p2 ' 

Rc = tl{l + 2pc cos J + pI+pI), 

Rn = il-tcf+tlipl+pl)-2tc{l~tc)pcC0S-f 



e(2t^p,)cos7+ lify/l-eH^ 2tc{l+PcC0Sj) , 



RnC^a^o = 2tcSin7 tcPsTPs\/'^-(^tl + epctc 



(6) 



The =F signs in the last two equations should be chosen 
to be the same, and correspond to whether the apex of 
the triangle in Fig. Q^) is to the right/left of the (0,0) 
point. Since both |Ac| and |A,j| are absorbed into the 
hadronic parameters Ps and pc, so there is no added un- 
certainty from \yii.h\- For the CKM angle /3 we use the 
latest average |ll[l2, /3 = 23.3° ± 1.5°. 

The full isospin analysis requires solving the five equa- 
tions © to obtain the parameters Pc, Ps, tc and e defined 
in Q and the weak angle 7. From St^+tt- and C^+^- 
one obtains two solutions for the parameters Ps and Pc 
as functions of the angle 7. Using these, Rc determines 




FIG. 2: Isospin analysis showing the hadronic parameters 
{pc,Ps,tc,t} versus 7 using current central values of the 
B — > TTTT data. Solutions for 7 occur at crossings of the 
curves. Experimental uncertainties are not shown, and are 
especially large for £3,4. This plot shows only one of two al- 
lowed (pcPa) solutions and one of the two allowed 7-regions. 



tcil)- Finally, i?„ and RnCT^o^^o each give two quadratic 
equations for e, which in general have four intersections 
in the e — 7 plane. We call ei.2 the two solutions from i?„ 
and 63,4 the two solutions from i?„C^o^o. An example of 
this GL isospin analysis is shown in Fig. |21 where we use 
the current central values for the data. For illustration 
we picked the solution for p^ and pg with \P/T\ < 1, but 
have shown all four e^'s. 

An obvious feature in Fig. [5] are the isospin bounds on 
7. It is well known that there are bounds on 7 in the 
absence of a measurement of C^o^o jl5l |. To find these 
analytically, one defines 7 = tt — /3 — acB + where 



sin(2a,ff) = 5„,(l-C72j-i/2 = _ 
cos(2^?)>(i?-l)(l-C2j-i/2 



-0.66 ±0.14, (7) 
= 0.53 ±0.19, 



with R = (l + i?c-i?„)^/(2i?c). The four solutions are 



163.° < 7 < -105.°, -31.i 
17.1° < 7 < 75.2°, 148. 



° < 7 < 26.3°, 
<7 < 206.°, (8) 



with uncertainty ±8.2° on each limit. At each of these 
bounds the two solutions ei^2 become degenerate and be- 
yond they are complex, indicating that the isospin trian- 
gle does not close. 

Solutions for 7 from the GL isospin analysis are given 
where the curves £3,4 intersect the curves ei_2- There are 
up to four solutions within each isospin bound, which 
are symmetric around 7cff = tt — /? — acff. We show in 
Fig. lithe results for 17.1° < 7 < 75.2°. The current 
central values £3^4 do not intersect ei.2, and in the ab- 
sence of experimental uncertainties there would be no 
solution for 7. The current central values for the ob- 
servables are such that the solutions for e from i?„ and 
i?„C^o^o are almost tangential. Including experimental 
uncertainties a large range of 7 is allowed, with the high- 
est confidence at 7 = 27° and 7 — 65°. This conclusion 
agrees with the CKMfitter group's analysis which incor- 
porates C^o^o [T^ . 

Using SCET at LO in as(mb) and AqcD/wh we have 
e = j5j], which corresponds to flat isospin triangles in 
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Neglecting EW-penguins, e is an RGE invariant quantity 
since Eq. ((HJ fixes it in terms of observables. Eq. Q 
makes an extraction of 7 from B — > tttt possible without 
needing precision data on In this method the 

central values for 7 are determined by finding the places 
where the ei and/or 62 curves cross the x-axis, meaning 
we solve £1,2(7) = 0- The other hadronic parameters, 
Pc Ps, and tc are determined in the same way as in the 
isospin analysis. This proposal for determining 7 using 
Eq. ^ is the main result of this letter. 

Using the central values for all the data besides C^o^o 
and solving £1.2(7) — gives the solutions 



7 --159°. 



-105° 



21.5°, 74.9°. 



(10) 



We have four solutions rather than the eight of the isospin 
analysis (which occur within the first and third isospin 
bounds in Q), because factorization for the _B — > tttt 
amplitudes resolves the discrete ambiguity in ps and pc 
in favor of |P/T| < 1 solutions (this follows from the 
factorization for light-quark penguins, the size of Wilson 
coefficients jCharm velocity power counting, and factors 
of as{mc) 0,0)- Next we analyze the theoretical and 
experimental uncertainties in our method for 7, and con- 
trast these with the isospin analysis, focusing on the two 
solutions which can occur in the 17.1° < 7 < 75.2° region 
preferred by global fits for the unitarity triangle [l^ . 
To estimate the theoretical uncertainty we take 



0.2 < £ < 0.2, 



(11) 



which corresponds to roughly a 20% effect from pertur- 
bative or power corrections. We also consider a much 
more pessimistic scenario where this range is doubled to 
£ = ±0.4. Note that |£| < 0.2 can accommodate the 
so-called "chirally enhanced" power corrections, which 
have been argued to dominate . Using the results from 
Ref. 0, including known as{mb) and power corrections, 
we randomly scan the two complex parameters Xa and 
Xh in the BBNS range to find |£| = |Im(C/T) |qcdf = 
—0.05 ± 0.04. This is below the uncertainty assigned to 
our analysis. 

In Fig. 131 we show £1,2 for the region 65° < 7 < 78°. 
Here the solution is 7 = 74.9°, and the different shading 
corresponds to the theory uncertainty with |£| < 0.2(0.4). 
The solution for 7 is very close to the isospin bound, so 
the upward uncertainty on 7 is very small. (The un- 
certainty in the isospin bound is contained in the ex- 
perimental uncertaintiy.) For the downward uncertainty 
we consider the overlap with the shaded region. For 
|£| < 0.2 we find A7theo =t1fo, while for |£| < 0.4 

we find A7thco =^5 2° ■ On top of that there are uncer- 
tainties from hadronic isospin violation, typically ^ 3%, 
which we take to be ±2° since no source of enhanced 
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FIG. 3: Regions of 7 preferred by our analysis. The shaded 
areas show our best estimate of the theoretical uncertainty 
from power corrections, —0.2 < e < 0.2 as well as the pes- 
simistic estimate —0.4 < e < 0.4. Experimental uncertainties 
are not shown. 

isospion violation has been identified. (A slightly larger 
~ 5° uncertainty was found in but using a smaller 
penguin amplitude. Larger isospin violation can be ac- 
counted for by scaling up the ±2° lower bound on our 
theory error.) Thus, with perfect data at the current 
central values we arrive at a theory uncertainty with 
|£| < 0.2 as A7thoo = ±2°. Repeating for the smaller so- 
lution at 7 = 21.5°, we find a larger theory uncertainty, 
A7thoo = ^440, since the £1,2 curves are flatter near this 
solution. 

To determine the 1-cr experimental errors, we use the 
program Minuit. Taking £ = and fitting to 7 and the 
four hadronic parameters we find 



7 = 2i.5°t^,:^: 



7 — ''i-y -10.6° 



(12) 



These uncertainties are purely experimental and are 
propagated with the assumption that the original input 
data are uncorrelated. If we instead set £ = 0.2 then we 



find 7 



73.3°Ii3';3„ and 7 
-0.2 gives 7 



30.7 



0-1-11.1° 



_7 20 J whereas fix- 
ing £ = -u:z gives 7 = 75.2°l^;f„ and 7 = 17.2°l^;^°o. 
Combining these numbers we obtain our final result for 
7 including all sources of uncertainty 



7 74.9° ± 2°I 



9.4° 

13.3° 



(13) 



Here the first error is theoretical, and the last errors are 
experimental where we picked the largest range obtained 
in varying £ = ±0.2. The theory error increases to A7 = 
-1 2° fo'' ^1^6 more pessimistic case e = ±0.4. The analog 
of (O for the lower solution is 7 = 21.5°t^;f„t"gr . 

The analysis presented here relies on the fact that a 
small value of |£| is allowed only for a narrow range of 
7. While this is certainly true given the current central 
values of the data, it is instructive to investigate how the 
quality of the analysis is affected if the data central val- 
ues change. For example, it could be that the value of £1 
never exceeds 0.2, increasing the uncertainties from the 
small-£ analysis significantly. A second extreme situation 
is that £ never reaches zero. To study these questions, 
we generate random sets of data using Gaussian distri- 
butions with the current central values and width of the 
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FIG. 4: Uncertainty in the isospin analysis from C^o^o (y- 
axis) vs. theoretical uncertainty from our new method (x- 
axis). The upper (lower) two plots use a ±0.5 (±0.25) uncer- 
tainty in £3,4. The plots on the left (right) correspond to the 
solution near the lower (upper) isospin bound. 

1-a uncertainties. We generate 10000 sets of "data" , and 
after imposing sin(2aoff) < 1 and cos(20) < 1 are left 
with 9688 sets. Of these 96% have solutions for e = 0. 
For ci we find 88% (70%) of the sets have the maximum 
value above 0.2 (0.4). It is only for these data sets that 
our analysis works. For £2 we find that 100% of the 
sets have their minimum below -0.4. Thus, the small e 
analysis works well in most cases. 

We can also study the uncertainty in our analysis, com- 
pared to the GL isospin analysis. Rather than perform- 
ing a full error analysis for the 9688 sets, we use the 
following approximation. We assume that experimental 



uncertainty ACt^o^^o dominates, and compare the result- 
ing uncertainty in the GL isospin analysis to the theo- 
retical uncertainty in our analysis, for cases where values 
of 7 exist with e > 0.2 as discussed above. The cur- 
rent AC^OttO = ±0.39 gives rise to a Ae^^i ~ ±0.5. In 
Fig. 21 we show the uncertainties in the GL analysis com- 
pared with the theoretical uncertainties of the analysis 
presented here, for both solutions of 7. The plots use 
4000 points. If we take Ae^^^ ~ ±0.25, the GL analysis 
still has uncertainties in 7 that are considerably larger 
than the small e analysis. We also see that the ±2° error 
quoted in 113|) is typical. 

We have presented a new method for obtaining 7 from 
B ^ TTTT decays without C^o^o . Our analysis uses SCET 
to eliminate one hadronic parameter. The theory uncer- 
tainty for a solution 7 = 74.9° are small, ±2 or I520J 
depending on the estimate for power corrections. Ana- 
lyzing possible future shifts in the data and decreases in 
the C^o^o uncertainty, we find that this method should 
have smaller uncertainty than the isospin analysis for 
quite some time. The analysis can be redone including 
the electroweak penguins. Obviously agreement between 
Babar and Belle on S'tt+tt-, C'tt+tt- is needed before one 
will have complete trust in the 7 from B tttt. 
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